THE MEAN-FIELD LIMIT FOR A REGULARIZED 
VLASOV-MAXWELL DYNAMICS 



FRANCOIS GOLSE 

Abstract. The present work establishes the mean-field limit of a 
A-particle system towards a regularized variant of the relativistic 
Vlasov-Maxwell system, following the work of Braun-Hepp [Comm. 
in Math. Phys. 56 (1977), 101-113] and Dobrushin [Func. Anal. 
Appl. 13 (1979), 115-123] for the Vlasov-Poisson system. The 
main ingredients in the analysis of this system are (a) a kinetic 
formulation of the Maxwell equations in terms of a distribution of 
electromagnetic potential in the momentum variable, (b) a regular- 
ization procedure for which an analogue of the total energy — i.e. 
the kinetic energy of the particles plus the energy of the electro- 
magnetic field — is conserved and (c) an analogue of Dobrushin's 
stability estimate for the Monge-Kantorovich-Rubinstein distance 
between two solutions of the regularized Vlasov-Poisson dynamics 
adapted to retarded potentials. 



1. Introduction 

Vlasov equations are kinetic equations used to describe the evolution 
of dilute, collisionless systems of particles coupled by some long-range 
interaction. For instance, in the case of a rarefied plasma, each charged 
particle is accelerated by the electrostatic field created by all the other 
particles, a situation described by the Vlasov-Poisson system [20]. A 
similar Vlasov-Poisson system is used in astrophysics to model the col- 
lective behavior of like massive objects coupled by Newton's gravita- 
tional force — which is attractive, at variance with the electrostatic 
force between charged particles with charges of the same sign. 

Other types of interactions can also be considered. In the case of 
magnetized plasmas - - e.g. plasmas in tokamak machines - - each 
charged particle is accelerated by the Lorentz force resulting from the 
electromagnetic field created by all the other particles, a situation de- 
scribed by the Vlasov-Maxwell system [13]. There are also general 
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relativistic variants of the gravitational Vlasov-Poisson system, such as 
the Vlasov-Einstein [5] or the (model) Vlasov- Nordstrom systems [4]. 

At the time of this writing, there is no rigorous derivation of any 
of these kinetic equations from the corresponding iV-body problems, 
in the large N limit. Braun and Hepp [3] and Dobrushin [7] have 
proposed instead rigorous derivations of a system analogous to the 
Vlasov-Poisson system, replacing the Coulomb potential with a twice 
differentiable mollification thereof. Hauray and Jabin [10] have suc- 
ceeded in treating the case of singular potentials, but their analysis 
does not include the Coulomb singularity yet. It should be noted that 
the situation is significantly better in the quantum case: the Hartree 
or Schrodinger-Poisson equation, that is the quantum analogue of the 
Vlasov-Poisson system has been derived from the quantum iV-body 
problem with Coulomb potential by Erdos and Yau [9]. 

In the case of the Vlasov-Maxwell system, it is not even clear what 
the corresponding iV-body problem should be. For instance, there is at 
present no satisfying description of the electromagnetic self-interaction 
(i.e. the action of the electromagnetic field created by a moving, point- 
like charge on itself) within the theory of classical electrodynamics. 
Bibliographic references discussing this well-known issue can be found 
in the remarks following the introduction of the regularized iV-body 
system (27). 

The present paper proposes an analogue of the Braun-Hepp or Do- 
brushin theory for (a mollified version of) the Vlasov-Maxwell system. 
Elskens, Kiessling and Ricci [8] have recently considered a special rel- 
ativistic variant of the gravitational Vlasov-Poisson system where the 
Poisson equation for the potential is replaced with a linear wave equa- 
tion. 

The problem of deriving a regularized variant of the Vlasov-Maxwell 
equations from a particle system is explicitly mentioned by Kiessling 
on p. Ill in his survey article [12]. As expected, this problem can 
be handled more or less along the line of [8], with however a few sig- 
nificant modifications, which it is the purpose of the present work to 
explain. As we shall see, the regularization procedure removes all the 
conceptual difficulties pertaining to the electromagnetic self-interaction 
mentioned above. The idea of considering a mollified variant of the 
Vlasov-Maxwell dynamics is therefore motivated by reasons other than 
analytical simplicity. 
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2. A SCALAR FORMULATION OF THE VLASOV- MAXWELL SYSTEM 

Regularizing the Coulomb potential is a crucial step in the mean- 
field limit established by Braun and Hepp [3] and Dobrushin [7]. At 
first sight, it is not clear that the same procedure can be applied to 
the Vlasov- Maxwell system, as the electromagnetic field involves both a 
scalar and vector potentials. Perhaps for that reason, Elskens, Kiessling 
and Ricci [8] chose "to purge the inhomogeneous wave equation for [the 
vector potential] A and all terms involving (derivatives of) A in the 
Lorentz force" and retain only the inhomogeneous wave equation for 
the scalar potential. 

In this section, we explain how both the scalar and vector potentials 
in the Vlasov-Maxwell system can be expressed in terms of a single 
scalar potential distributed in the momentum variable. This formu- 
lation, introduced in [1, 2] for other purposes — viz. for studying 
regularity properties of the Vlasov-Maxwell system — although per- 
haps not absolutely necessary in the present context, greatly simplifies 
our arguments in the sequel. 

The Vlasov-Maxwell system for a single species of particles with 
momentum distribution function f(t,x,C,) — meaning that f(t,x,£) is 
the density of particles with momentum (eR 3 that are located at the 
position x G R 3 at time t > — is written as follows: 



dtf + v(0 ■ V x f + (E + v(0 A B) ■ V^/ = , 

div x B = , 

d t B + curl s E = 0, 

div x E = p f , 

d t E - curlz B = -j f . 



We have used the notation 



v(0 ■= V € e(0 = J— , where e(£) := ^lTW 
V 1 + If I 



and 



p f (t,x):= [ f(t,x,Z)<%, while j f (t, x) := / v^)f(t,x,^. 

</R3 JR3 

This system is posed in the whole Euclidian space — i.e. for all 
x, £ G R 3 — and supplemented with the initial data 

(2) /U = / m > E Lo = Em > 4=o = ' 
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where E m is the electrostatic field created by the distribution of charges 

Jin. 



(3) 



E % 



-V x <j> m , <p m 



R3 



(We restrict our attention to such initial data for simplicity; treating 
the case of general initial data 

/U = r> E Lo = E ' m > B Lo = B ' m > 

with the compatibility conditions 



div x E m = / f m d£, div x B m = 

is not significantly more complicated, as it only involves one extra 
linear, homogeneous wave equation.) 

Solve for <po the Cauchy problem for the linear wave equation 



( n ttS< i> = o , 



(4) 



< 0o| t=o = m = (-A :c )- 1 



R3 



I <9 t 0oU = O. 

The scalar formulation of (1) is as follows. Consider the coupled 
system with unknowns / = f(t, x,£) > and u = u(t, x,£) G R: 

f dtf + v(Z)-V x f + (E + v(t)AB)-Vzf = 0, 



(5) 



= 00+ / udi, A= v(0ud£, 

JR3 JR3 



[ B = curl, A , E = —d t A - V x , 
with initial data 

( 6 ) /U = / iB . «U = ft«U = «- 

Physically, u(t,x, £) is the instantaneous Lienard-Wiechert potential 
created by a particle with momentum £ at the position distributed 
under f(t,x,g). 

Proposition 2.1. Let / e C c °°([0,T] x R 3 x R 3 ), and let O be the 

solution of (4)- If (f,u) satisfies (5) with initial data (6), then the 
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electromagnetic potential (0, A) defined as in (5) is smooth and satisfies 
the Lorentz gauge 



while (f,E,B), with the electromagnetic field (E,B) as in (5), satisfies 
(1) with initial condition (2)-(3). 

We shall use systematically the following elements of notation: C C (X) 
(resp. C*(X), C%°(X)) designates the set of continuous (resp. C k , C°°) 
functions with compact support included in X. Likewise, Cj,(X) des- 
ignates the set of bounded continuous functions defined on X, while 
C*(X) designates the set of C k functions defined on X whose deriva- 
tives of order < k are bouned on X. Finally, C£°(X) designates the 
class of C°° defined on X whose derivatives of all orders are bounded 
on X. 

Proof. Since f m G C£°(R 3 x R 3 ), the initial electrostatic potential 



so that the solution O of the Cauchy problem (4) is in C£°([0, T] x R 3 ) 
for each T > 0. 

Since / G C£°([0, T] x R 3 x R 3 ), the finite speed of propagation and 
the regularity theory for solutions of the wave equation implies that 
u G C c °°([0,T] x R 3 x R 3 ). 

Therefore the scalar potential G C£°([0,T] x R 3 ) and the vector 
potential A G C£°([0, T] x R 3 ; R 3 ), so that, in view of the last equalities 
in (5), the electric field E G C fe °°([0,T] x R 3 ; R 3 ) while the magnetic 
field B G C c °°([0, T] x R 3 ; R 3 ) and 



Besides, p f and j f belong to C c °°([0, T] x R 3 ) and C c °°([0, T] x R 3 ; R 3 ) 
respectively, and since v(£) = V^e(^), one has div^(v (^) A B{t,x)) = 
B{t, x) ■ curl € v(£) = so that 



d t (j) + div x A = 



<F = (- 



A^- 1 ^ e C 6 °°(R 3 ) , 



div x 5 = 0, 

d t B + cur\ x E = . 
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Let us verify that (0, A) satisfies the Lorentz gauge condition. Setting 
g := d t <j> + d\v x A G C b °°([0, T] x R 3 ), one has 

U t>x g = D t)X d t (p + U t!X / (d t u + dxv x (v(£)u))d£ 

= d t p f + div x j f = 0, 
together with the initial conditions 

g\ t=0 = ^t0o| t=o + / (d t u\ + div x (v(£)u\ ))d£ = , 

</R3 

and 

d t g\ t=0 = d^\ t=0 + div x d t A\ t=0 

= n tjX (/) \ +A x (f)\ + / d 2 t u\ d£+ I v((,)-V x d t u\ d£ 

JR 3 iR3 



= A x <j> m + [ n tiX u\ t=0 d£+ [ A x u\ t=0 dt 

J R 3 ./R 3 



= A x m + / f"^ = 0. 
</r 3 

By the uniqueness property for the wave equation, one concludes that 
g = 0, which is the Lorentz gauge condition. 

It remains to verify that the electromagnetic field satisfies the Gauss 
and Maxwell-Ampere equations: by the Lorentz gauge condition, 

div x E = D tjX (f> - d t (d t (j) + div x A) 

= D t>x ^ + / ud£= fd£ = p f , 

JR3 JR3 

while 

d t E - curl, B = -Dt >x A - V x {d t 4> + div x A) 

= -U t>x [ v(0ud{; = - [ v(t)fd£ = -J f • 

JR3 JR 3 

Therefore (/, E, B) is a solution of the Vlasov-Maxwell system (1), and 

E Lo = -[ <mu\ t=0 di-v x ! u \ t=0 dt-v x <f> \ t=0 

</R3 JR 3 



= -V x (-A x y 1 / f m dt = E m , 

</R3 



while 

B 



= curL J ^ v(^)u\ t=0 d^ = , 
so that (E,B) also satisfies the initial condition (2). □ 
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With this formulation of the Vlasov-Maxwell system, the analogue 
of the Coulomb potential in the Vlasov-Poisson system becomes fairly 
obvious. Let Y be the forward fundamental solution of the d'Alembert 
operator, i.e. the only Y G V'(R x R 3 ) satisfying 



(7) 



O t , x Y = 8(t, x )=(o,o) , 
supp(V) c R+ x R 3 



We recall that, in space dimension 3, the distribution Y is given by the 
formula 

(8) Y=^5(\x\-t), 

where the notation S(\x\ — t) is understood as the surface measure on 
the sphere of radius t > centered at the origin. Then, the function u 
in (5) is given in terms of / by the formula 

(9) <,.,0 = ^* t ,*(W(-,-,0), £eR 3 

where * tjX denotes convolution in the variables (t,x), while the scalar 
potential O is given in terms of f m by 

(10) Mt,-) = dtY(t,-)* x G* x [ f in (;0dC, t>0, 
where 

(11) G{z) = -±- zeR 5 

4:7l\Z\ 

is the fundamental solution of —A x . 

Therefore, the Lorentz force field F := E + v A B in the Vlasov- 
Maxwell system is given in terms of the particle distribution function 
/ by the formula 

F[f] = - I v x d t Y{t,-)* x G* x f$,;-n)d-n 

</R 3 

/ {^x + v(v)dt)Y* t>x f(;;7l)dri 

" / v{0 A (v(rj) A V X Y M, x f(; ; r,))dr] , 

</R3 

to be compared with the formula 

i^PoissonL/] = - / V X G * x f(t,-, f])df] 
^R3 

in the case of the Vlasov-Poisson system. 
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In the Vlasov-Poisson case, the iV-particle dynamics considered by 
Hauray and Jabin (see for instance formula (1.1)-(1.2) in [10]) is the 
differential system 



(12) xtf) = -jfzr~[ V - G (^( t ) - X M , 1 < i < N ■ 

3 = 1 
3& 



By analogy, in the Vlasov-Maxwell case, one should consider the 
system 
(13) 

{ ±i(t) =v(£i(t)) 

1 N 

= ~w^i E V *W> •) *x C(^(0) - ^(0)) 



3 = 1 

37^ 



N —Y, J Q ^ x + v^(s))d t )Y(t~s,x l (t)-x 3 (s))ds 

i N r l 

3 = 1 ^ 



i = l,...,iV. 



There is an obvious difficulty in (12): since V X G is singular at the 
origin, the right-hand side of the equation for xi is not defined whenever 
Xi(t) = Xj(t) for some j ^ i. More seriously, in (13), the second and 
third terms in the r.h.s. of the equation for £j are not a priori well 
defined quantities because Y is a measure and not a (C 1 ) function. Yet 
this difficulty disappears in the reference [3], where the Green function 
G is replaced with a regularization thereof. Likewise, we shall replace 
Y with some regularized variant thereof, which removes the difficulty 
mentioned above. 

Notice that, while (12) is a system of N coupled ordinary differential 
equations, (13) is a system of N coupled integro-differential equations 
- more specifically, delay differential equations. The same is true 
of the dynamical system (7)-(9) considered in [8], and this is a conse- 
quence of the retarded potential formula for the solution of the Maxwell 
system, or, equivalently, of Kirchhoff 's formula (9) for the solution of 
the wave equation. 
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3. A REGULARIZATION OF THE VLASOV-MAXWELL SYSTEM 

Any regularization of Y would transform (13) into a well-posed sys- 
tem. Yet, classical solutions of the original Vlasov-Maxwell system 
satisfy some conservation laws, which one might wish to preserve by 
the regularization procedure. 

Since v = V^e, 

(14) div Xi £(u(f ), E + v A B) = div^v A B) = B ■ curl ? v = 

so that the characteristic field of the transport equation in (1) preserves 
the phase space Lebesgue measure dxd^. Hence, if (/, E, B) is a classi- 
cal solution of (1) on the time interval [0,T], then, for all p G [1, +oo], 
one has 

(15) •, -)|Uf(R 3 xR3) = ||/ m (-, -)IUf(R 3 xR3) 

for each t G [0,T]. This obviously remains true if the electromagnetic 
field (E, B) is replaced in the transport equation governing / with any 
regularization thereof in the t, x variables. 

Moreover, if (/, E, B) is a classical solution of (1) on the time interval 
[0,T] with, for instance, / G C£([0,T] x R 3 x R 3 ) while the electro- 
magnetic field E,B G Cl([0,T] x R 3 ; R 3 ) fl C([0, T]; L 2 (R 3 ; R 3 )), then 

(16) II e(£)f(t,x,£)dxd£ + \ f (\E\ 2 + \B\ 2 )(t,x)dx = Const., 

which is the conservation of energy (the first integral being the kinetic 
energy of the particles while the second is the electromagnetic energy). 
Conservation of energy is not expected to be preserved by all regular- 
izations of the electromagnetic field. 

There is however one clever regularization procedure for the Vlasov- 
Maxwell system that preserves a quantity analogous to the energy and 
which converges to the energy (16) as the regularization is removed. 
This regularization of the Vlasov-Maxwell system, proposed by Rein in 
[16], is recalled below for the reader's convenience 1 . 

Let x e Q°(R 3 ) satisfy 

X (x) = X (-x) > , supp(x) C 5(0, 1) , / X (x)dx = 1 , 

JR3 

and define the regularizing sequence 

Xe(x) = ^X (^) 



Prof. Rein kindly informed the author that this regularization procedure had 
been introduced earlier in this context by E. Horst in his Habilitationsschrift. 
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for each e > 0. Consider then, for each e > 0, the following regularized 
variant of (5) with unknown (f e ,u e ): 

{ d t f e + v(0 ■ V x f e + (E e + v(0 A B e ) • V € / e = , 

I— \,x^e Xe *x Xe *x fe 5 



(17) 



<Pe = Xe *x Xe *x 00 + / U e d£ , A e = / v{£)ll e d£ , 

'R 3 iR 3 



, B t = cml x A t , E t = -d t A t - V x <f) t , 
where 0o is defined in as (4), and with the same initial data as for (5): 



(18) 



fe 



t=0 



t=0 



d t u e 



0. 



Along with the electromagnetic field (E e: B e ) of the regularized sys- 
tem (17), consider (E e , B € ) defined as the solution of the Maxwell sys- 
tem 



div^ B e = 



(19) 

with initial data 
(20) E e 



d t B t + curl, E t = , 

div x E € =Xe*xPf e , 

d t E e - curU B e = -Xe * x jf e , 



t=0 



Xe*xE m , Be\ t=Q = 0. 



Proposition 3.1 (Rein [16]). Let e > and f in G L°°(R 3 x R 3 ) be 
such that 

f m > a.e. on R 3 x R 3 and supp(/ m ) is compact. 

Then the regularized system (17) has a unique weak solution (f e ,u e ) 
defined for all positive times and satisfying u t G C%°([0,T] x R 3 x R 3 ) 
for each T > 0, while 

fe(t,x,£)>0 for a.e. (t, x, f) G R+ x R 3 x R 3 , 
||/ e (i, •, -)IUf(R 3 xR 3 ) = Const, for all p G [l,+oo], and 



IL 



'R s xR» 

for each t > 0. 



e(Z)f e (t,x,Z)dxd£ + ± / ( | E e | 2 + \B t \ 2 )(t, x)dx = Const. 



R 3 
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Proof. The same computation as in the proof of Proposition 2.1 shows 
that (/ e , E e , B e ) satisfies 



(21) 



{ dtfe + v(0 ■ V x / e + (E e + v(0 A S e ) • VJ e = , 
div x B e = , 
<m + curl, £ e = , 

dW x E e =Xe *x Xe *xPf e , 
d t E € - CUl\ x B e = -Xe *x Xe *x 3 f € , 



with initial data 

(22) fe\ t=Q = f m , E e \ t=Q = Xe*xXe*xE m , B t \^ = , 

where E m is defined as in (3). 

This is precisely the regularized Vlasov-Maxwell system studied by 
Rein in [16], to which we refer for the existence of a solution. 

That f e > a.e. on R + x R 3 x R 3 and t \\f € (t, •, •) ||lp(r 3 xr») 
is constant on R + follows from the method of characteristics for the 
transport equation satisfied by f e since E e and B e are smooth with 
bounded derivatives of all orders on R + x R 3 x R 3 and the vector field 

(x, H> ^(i, x) + u(0 A B e (t, x)) 

is divergence free, while f m > a.e. on R 3 x R 3 and f m G L P (R 3 x R 3 ) 
for each p £ [1, +oo]. 

We briefly recall from [16] the elegant argument leading to the con- 
servation of the (pseudo-) energy for (21), which is the reason for the 
specific regularization procedure chosen in (21). For simplicity, we as- 
sume that f m G C^(R 3 x R 3 ), so that, by the method of characteristics, 
f e G Cl([0,T] x R 3 x R 3 ) for each T > 0. Then 

(23) ±[[ e(t)f e (t,x,t)dxdt= f E e -j fe (t,x)dx. 

By uniqueness of the solution of the Maxwell system, E e = Xe*x E e , so 
that 

/ E e ■ j ft (t, x)dx = / (xe * x E e ) ■ j fe (t, x)dx 
= E e - (xe* x jfj(t,x)dx 

JR3 

where the last equality follows from the fact that Xe is even. On the 
other hand, the classical computation leading to the conservation of 



(24) 
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energy in the Maxwell system (19) is as follows: 

^ / \{\E e \ 2 + \B e \ 2 ){t,x)dx 
at 7 R 3 

= / (E e -d t E e + B e -d t B e )(t,x)dx 
= - / E e - ( X e*x3f e )(t,x)dx 
+ / (E e ■ cur\ x B e — B e ■ cur\ x E e )(t, x)dx . 

J-R3 

Now the last integral above vanishes since 

E t ■ cur\ x B e — B e - cur\ x E e — — div x (E t A B e ) 

and B t (t, •) has compact support for each t > 0. Combining (23) and 

(24) give the announced conservation of (pseudo-) energy for (21), and 
therefore for (17). □ 

At this point, we introduce the regularized iV-particle dynamics for 
the Vlasov-Maxwell system. Define 

(25) Y e = X e*xXe*xY. 

In other words, Y € is the solution of 

O t , x Y e = 5 t=0 ® ( X e * x Xe) , in V(R x R 3 ) 
supp(F e ) C R+ x R 3 

or, equivalently, ^e| R+xR3 is the solution of the Cauchy problem 

n ttX Y e = , x e R 3 , t > , 

d tY e \ t=0 =Xe *x Xe- 

The classical theory of regularity and the finite speed of propagation 
for solutions of the wave equation imply that, for each e > 0, 

(26) Y e G C°°(R + x R 3 ) and supp(F e (t, •)) C 5(0, t + 2e). 
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Next we consider the system of (delay) differential equations 
(27) 

' ±i(t) = v(£i(t)) 
1 N 

i=i 
f * /•* 

< ~mT, (^ x + v(^(s))d t )Y e (t-s,x i (t)-x j (s))ds 
1 * /•* 

~mJ2 "feW) A W^))Av I ^-^,w-i 3 ( S ))) ( i s . 

i = l,...,iV. 



Like (13), this is a system of N coupled delay differential equations. 
However, unlike in the case of (13), the r.h.s. of the second equation 
above is obviously well-defined since Y e e C£°([0, T] x R 3 ) for each 
T > 0. 

Comparing the original dynamics (13) with its regularized analogue 
(27), one notices that, in (13) the force acting on the i-ih particle 
is the electric or electromagnetic force created by the N — 1 other 
particles. Self-interaction — i.e. the action of the electromagnetic field 
created by the i-th particle on itself — is therefore neglected in (13), 
unlike in the regularized model (27). As a matter of fact, there are 
serious conceptual difficulties with the notion of self-force in classical 
electrodynamics, which are beyond the scope of the present paper. The 
interested reader is referred to the discussion on pp. 675-676 of [8] - 
see also chapter 16 of [11] and especially [17] for more on this subject. 
That the self-interaction force is present in the regularized model (27) 
is only for mathematical convenience (see below). 

In fact, if the regularizing parameter is kept fixed, neglecting the self- 
interaction force in (27) will only produce an error of order 0(1/N) over 
finite time intervals. 
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Indeed, consider instead of (27) the following regularized variant of 
(13): 
(28) 

( Xi(t) = v(ii(t)) 

1 N 

= -jyTi E ^dtY e (t, •) *x G(£i(0) - xj(0)) 



3 = 1 
3& 



i " r* 

-J^jY, J o (^ X + v(Cj(s))d t )Y e (t - s, Xi (t) - xj(s))ds 

1 N /•* 

--^ / 4(t))A(4-W)AV^ e (t-Mi(t)-i;(«)))<is. 

3 = 1 ^° 



i = l,...,JV. 



That the regularized dynamics defined by (27) and (28) are equiva- 
lent on finite time intervals as N — > +oo is summarized in the following 
statement. 



Proposition 3.2. Let e, T > and N > 2. If (x i: £i)i<i<N and 
( x i,£i)i<i<N are solutions of (27) and of (28) respectively on [0,T], 
and if 



Xi (0) = Xi(0) , £(0) = 6(0) fori = l,...,N, 



then one has the error estimate 



\xi(t) - Xi (t)\ + \ii(t) - &(t)\ = 0(1/N) 



for alii = 1, . . . , N and all t E [0, T] . 



MEAN-FIELD LIMIT FOR VLASOV-MAXWELL 

Proof. Observe that, for each N > 2 and 



15 



i N r* 

M —V] / (v x + u(4(s))^)y e (*- 3,^(0-^(3))^ 

iV ~ 1 j=i J o 

i N r* 

7V 3=1 J ° 
1 N /•* 



- (V, + u(&(s))$)y e (i - s, ^(i) - ^(s))|ds 

1 /■* 



2£ 

< jy _ ]_ ll^t,a:^IU°°([0,r]xR 3 ) 
JV 

! - 

+ 



1 /" 

^E H V M y e||L~([0,T]xR3) ^ l%( S ) -^(s)|ds 
1 " _ /■* 

+ 2T||v^y|| L ^ ([0 ,T]xR3)l^(t) - • 



The difference between the right-hand sides of the equations for & 
in (27) and (28) involves two more terms analogous to this one. 
Thus, denoting 



Si(t) :=\Ut)-Xi(t)\ + \£i(t)-Zi(t)\, 
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we see that, for each i = 1, . . . , N and all t G [0, T], one has 

(ft i N ft fT a. _|_ 9 

/ S i (r)dT + -J2 / SitfdsdT + 

= C e / 5i(r)dr + -X; / (f - a)£(a)da + 
/ T N rt 



N 
4T + 2 



TV 

since £j(0) = 0, where C e is a constant depending on 

||Vt jX l^||£<x>([o,T]xR 3 ) ) II Vi j2 ,F e ||L°o([o,T]xR 3 ) , ||C*s V x l^ ||l°o([o,T] x R 3 ) 

and ||V^|| L <x, (R 3) . 

Therefore, summing for % = 1, . . . , N both sides of the inequality 
above shows that 

1 N ( /"* 4T -l- 9 



N 

and we conclude from Gronwall's inequality that 

S{t) < C e ^le c ^ 1+T » , < t < T . 

N ~ ~ 



Substituting this in the inequality satisfied by Si, we find that 

4T + 2 S 
N 

2 
N 

so that, applying again Gronwall's inequality leads to 



Si(t) < C e ^jf Si(T)dr + T £(s)ds + 
<C e jf ^(r)dr + — e 



C e (l+T)t + 



< -^(1 + 2Te Ce(1+T)T )e at , 0<f<7\i = l,...,JV. 



□ 



4. Dobrushin's estimate 
Let A; be a Revalued Radon measure on R 1+d x R 1+d of the form 
(29) k(t, z, r, C) = m(t, z, ()5 T=0 + r(t, z,r,C), z, ( G K d , s, t > . 
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We assume that r G C b (R 1+d x R 1+d ; R n ) and m G C b (R 1+d x R d ; R n ) 
satisfy the properties 

r(t, z, r, () = whenever r [0, t] and z, ( G K d , 

m(t, z, C) = whenever t < and 2,(6^, 

together with the Lipschitz condition 

z, r, C) - r(t, z', r, (') \ + \m(t, z, Co) - r(t, z', Q \ 

< Up T (k)(l A \z - A + 1 A |C - C'l + 1 A |Co - Col) 

for all t, t 1 ', r, r' G [0, T] and all z, z', C, C'> Co, Co e R d > where the notation 
a Ab designates min(a, b) for all a, b G R. 

We consider the integral operator with kernel k, i.e. 

Kfi(t,z)= I r(t,z,T,()n(r,d()dT+ I m(t, z, C)//(0, dC) 

for each /i G C(R+; u> — "P(R d )), where w — "P(R d ) designates the set 
of Borel probability measures on R d equipped with its weak topology. 

For any given Borel probability measure p m on R d , consider Z = 
Z(t,z;p m ), the solution of the mean-field integro-differential equation 

r z = K P (t,z), 

(32) I p(t,-) = Z(t,-;p m )#p m , 

( Z(0,z;p in ) = z. 

(If $ : R d i — y R d is a Borel transformation and p a Borel probability 
measure on R d , the notation $#p designates the image measure of p 
under <3>, defined by the formula $#p(A) = p($ _1 (A)) for each Borel 
set A c R d .) 

Proposition 4.1. For eac/i Borel probability measure p m on H d , the 
mean-field integro-differential equation (32) has a unique global solution 
Z G C(R + xR d ;R d ) such thatt i-> Z(t,z) is continuously differentiate 
on R + for each z G H d . 

Proof. To any bounded map Z belonging to C(R+ x TV 1 ; R d ), we asso- 
ciate the map Z in C(R + x H d ; H d ) defined by 

Z(t, z) = z+ f K(Z(; -)#p m )(r, Z(r, z))dr 
Jo 

= z+f f f ris,Z( S ,z),T,Z(TX))p m (d()dTds 
Jo Jo Jn d 

+ f [ m(s,Z(s,z),()p m (dOds. 
Jo Jn d 
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Let Z U Z 2 e C(R+ x K d ;K d ) and denote Z U Z 2 the maps so denned 
in terms of Z\, Z 2 . Then, by the Lipschitz condition (31) 

\Z x (t,z)-Z 2 (t,z)\ 

< Lip T (&) f f f \Zi(s,z) - Z 2 (s,z)\p in (d()dTds 
Jo Jo Jn d 

+ Up T (k) f f [ iZ^O-Z^OiridQdrds 

JO JO J~R d 

+ Lip T (/c) f [ \Z l (s,z)-Z 2 (s,z)\p m (d()ds, 
Jo Jn d 

for all t e [0,T]. Denote 

£(t) := sup \Zi(t,z) - Z 2 (t,z)\ , 
zeR N 

S(t) := sup \Z x {t,z) - Z 2 (t,z)\ , 
zeR N 

for all t > 0. Then, for all * G [0, T], one has 

£(t) < Lip T {k) ^jf (1 + s)S{s)ds + £{r)drds^j 

(33) = Lip T (A;) ^jf (1 + s)£(s)ds + (t - T)E{r)d^j 

= Lip T (ife)(l + t) / S(s)ds. 
Jo 

With this estimate, we construct the solution by the usual Picard 
iteration procedure. Define recursively a sequence of maps (Z n ) n > by 

Z n+ \t,z) :=z+ f K(Z n (;-)#p m )(T,Z n (r,z))dT, t>0, 
Jo 

Z°(t,z) :=z. 
Applying (33) iteratively with the notation 

£ n (t) = sup \Z n+1 (t,z) -Z n {t,z)\ 
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shows that, for all t G [0,T], one has 



(34) 



S n {t) < Lip T (fc)(l +t) [ S^it^dh 

Jo 

< Lip T (fc) 2 (l +t) I + I 1 £ n - 2 (t 2 )dt 2 dt 1 

Jo Jo 

< Lip T (A;) 2 (l + t) 2 [ ^ e n - 2 (t 2 )dt 2 dt x 

Jo Jo 



<Lip T (k) n (l + t) n - sup £°(t n ) 
nl o<t n <t 



since 



/ 



Since 



lo<t„<t n - 1 <...<t 1 <tdt n dt n _i . . . dti — — - 

n\ 



Lip T (A;) n (l +t) n t n 
> , < +oo 

n>0 



for each t > 0, we conclude that the sequence Z n converges uniformly 
on [0, T] x H d for each T > towards a solution of the integral equation 

(35) Z(t,z) = z+ f K(Z(;-)#p in )(T,Z(r,z))dT 

Jo 

defined for each t > and such that Z G C(R+ x R d , R d ). 

If Zi and Z 2 are two solutions of that integral equation, denoting as 
above 

£(t) = sup \Z x {t,z) - Z 2 (t,z)\ , 
zeR d 

we conclude from (33) that 

£(t) < Lip T (ifc)(l + i) / S(s)ds, 

Jo 

and, following the argument in (34), that 

sup S(t) < Lip T (A;) n (l + T) nl — sup £{t) . 

0<t<T n\ 0<t<T 

Letting n — > +oo while keeping T > fixed shows that £{t) — for 
each t G [0, T]; and since T > is arbitrary, Zi = Z 2 on R + x R '. This 
establishes the uniqueness of the solution Z of the integral equation 
(35). 
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Finally, since the integral kernels r and m are bounded continu- 
ous functions of all their arguments while p m is a probability mea- 
sure on R d , applying the dominated convergence theorem shows that 
t H> K(Z(-, -)#p in )(t, z) is continuous for all z G R d , so that t H> Z(t, z) 
is continuously differentiable in t for all z G R d and satisfies the differ- 
ential form (32) of the integral equation (35). □ 

Once the existence and uniqueness of the solution of the mean-field 
equation (32) is known, the next natural question is whether its depen- 
dence in the initial probability measure p m is continuous. The right 
topology in which to measure the proximity of two initial probability 
densities p % \-,p % 2 is obviously that of the weak convergence of probabil- 
ity measures on R d . An extremely convenient way to obtain quanti- 
tative information in that topology is to use the Monge-Kantorovich- 
Rubinstein distance, whose definition we recall below. 

Definition 4.2. Let V(R d ) be the set of Borel probability measures on 
R d . The Monge-Kantorovich-Rubinstein distance on V(R d ) is defined 
by the formula 

dist MKR (p, v) = inf / / 1 A \x - y\n(dxdy) 

where H(p, v) is the set of Borel probability measures on R d x R d whose 
marginals are p and v, i. e. 



IL 



(j)(x)i/j(y)7r(dxdy) = / <p(x)p{dx) / i/j(y)u(dy) 

R d xR d JR d Jn d 

for all bounded continuous functions (f>,ip defined on H d . 

This distance is also called the Wasserstein distance — although it 
was considered more than ten years earlier by Kantorovich and Rubin- 
stein, and appears in the work of Monge. 

Our (slight) generalization of the Dobrushin estimate presented in 
[7] is stated in the following proposition. 

Proposition 4.3. Let p™ and p™ be two Borel probability measures on 
H d , and let Zi, Z2 be the corresponding solutions of (32), denoted 

Zj{t, z) := Z(t, z; pf) for j = 1,2, t> 0, and z G R d , 

whose existence and uniqueness is established by Proposition 4-1, and 

Pj (V):=^(V)#Pr> for j = 1,2. 
Then, for each t > 0, 

dist MKR (pi(t, -),P2(t, •)) < (1 + tLip,(A;)) e ' 2Lip ^ distMKiM"^") • 
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The main difference between this result and Proposition 4 in [7] is 
that the proposition above applies to an integro-differential equation 
of the form (32) satisfying the causality assumption (30), whereas the 
differential system considered in Dobrushin's original contribution [7] 
did not involve memory effects. This is only natural since Dobrushin 
was interested in the Vlasov-Poisson equation, a non-relativistic model 
in which the force field is created by the instantaneous distribution of 
charges. 

There are other minor differences: for instance, unlike Dobrushin's 
the estimate obtained here is local in time and uses the Gronwall in- 
equality; besides, the integral operator K in (32) considered in the 
proposition above is not necessarily a convolution operator as in [7]. 
But the argument is mostly the same as in [7], and we give it below 
only for the sake of being complete. 



for each t > and Zi,z 2 G R d . In the second equality above, we have 
used the identity 



Proof. One has 



K Pl (t, Z ± (t, Zl )) - Kp 2 (t, Z 2 (t, z 2 )) 




(36) f <KZ j (T,Q)p?(dQ= [ 0(Q) Pj (T.dQ 
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for all bounded <j>, Borel measurable on R d , since Pjij, •) = Zj(t, -^pj 1 
— see the definition of the image of a probability measure under a Borel 
map following (32), that is equivalent to (36) with <p = 1 A . 

Pick 7r m to be any probability measure on R d x H, d with p™ and p l 2 
as marginals: one has, for each t > and z±,z 2 G H d 



K Pl (t, Z^t, Zl ))-Kp 2 (t, Z 2 (t, z 2 )) 



[ [[ (r(t,Z 1 (t,z 1 ),r,Z 1 (r,Ci)) 
'o J Jn d xR d 



- r(t, Z 2 (t, z 2 ),r, Z 2 (t, C 2 )))7r m (dCi, <K 2 ) 
+ [[ (m(t, Z ± (t, Zl ), Ci) - m(t, Z 2 (t, z 2 ), C 2 ))vr m (rfCi, ^2) ■ 

Hence, as a consequence of (31), for each t G [0,T] and z±,z 2 G R d 

\K Pl (t, Z 1 (t,z 1 ))-Kp 2 (t,Z 2 (t,z 2 ))\ 

<Lip T (k) (tl A\Z 1 (t,z 1 ) - Z 2 (t, z 2 )\ + V(r)dT + V(0^j , 

with the notation 

V(r):= [[ 1 A I^^CO-^r, (2)1^(^1,^2)- 

J </R d xR d 

Therefore, for each t G [0, T] and each z 1: z 2 G R d , one has 
|Zi(Mi) - Z 2 (t,z 2 )\ 

<\zi-z 2 \+ / \Kpi{s, Zi{s,zi)) - Kp 2 (s, Z 2 (s,z 2 ))\ds 
Jo 

< W ~ z i\ + Lip T (£;)t£>(0) 

+ Lip T (/c) ^jf si A \Zi(s, Zi) — Z 2 (s, z 2 )\ds+J^ V^drd^j . 
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Integrating with the measure n m (dzidz 2 ) both sides of the resulting 
inequality for 1 A \Zi(t, z±) — Z 2 (t, z 2 )\, we arrive at 

V(t) < (1 + tLip T (k))V(0) 

+ Lip T (k) ( [ sV(s)ds + [ [ V(r)drds 
\Jo Jo Jo 

= (1 + tUp T (k))V(0) 
+ Lip T (k) sV(s)ds + J (t - s)V(s)d 

= (1 + tLip T (k))V(0) + Lip T (A;)t f V(s)ds . 

Jo 

This implies the estimate 

V(t) < V(0)V(t) , t>0, 

where 

tt(f) = (1 + tLip t (A;))/ Lip < (fc) , t > . 

Now 

V(t)= [[ lA|Ci-C2k(t;dCidC 2 ) 

where 

n(t,-):=(Z l (t,-),Z 2 (t,-))#7r m 

is the image under the map (z±, z 2 ) (Zi(t, z±), Z 2 (t, z 2 )) of the mea- 
sure 7i m . Obviously 

n m e ii(p?,p 2 n ) =► tt(v) g n( Pl (t,-),p 2 (t,-)) 

for each t > 0, so that 

dist M ^(pi(V),P 2 (V)) < t>(t)T>(0). 

Minimizing the right-hand side of this inequality as n m runs through 
n(pl",p 2 n ) leads to the desired estimate. □ 



5. Application to the regularized Vlasov-Maxwell 

dynamics 

Let e > be a fixed regularization parameter, and consider the 
regularized Vlasov-Maxwell system (17). Henceforth we assume that 
f m is a probability density on R 3 x R 3 . Applying the method of 
characteristics to the transport equation governing f e shows that 

(37) f e (t, X e (t, x, v- f n ), E e (t, x, v- f m )) = f in (x, v) 
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where t >->■ (X e , E t )(t, •, •; f m ) is the solution of 

(38) I E e = E e (t,X e )+v(~ e )AB e (t,X e ) 

{ =F e [f e (t,;-)}(X e ,E e ), 

with the notation 

F e \f]:=- [ d t V x Y e (t,-)* x G* x f(0,; V )dr) 

(39) -/ (V x + ^)d t )Y e * tiX /(-,-, 77)^77 

Since div^^O, # e (t, x) + A B e (t,x)) = 0, the flow (X e ,E e ) 
leaves the Lebesgue measure dxdS, invariant so that (37) can be recast 
as 

(40) f € (t,-,-)dxdZ = (X e ,~ e )(t,-,-, f m )#f m dxdt . 

According to (39), the vector field (v(£), E e (t, x) + v(£) A B e (t,x)) 
can be expressed by an integral operator with kernel of the form (29) 
satisfying the causality and Lipschitz conditions (30) and (31): 
(41) 

(v(0, E e (t, x) + w (0 A B e (t, x)) = K e ((X e , E e )(t, •, •; /*»)#/*»)(*, x, £) 
= / // r e (t,x,£,T,y,r})f e (T,y,r})dydr}dT 

Jo JJrSxR 3 

+ / / m e (t, x, £, y, v)f in {y, v)dydr] . 

J iR 3 xR 3 

In this formula, r e = (r*,rf) and m e = (m\,m 2 e ), where r\, rf , m\ and 
ml are given by 

r\(t,x,£,s,y,rj) = 0, 

ml(t,x,£,y,ri) = lt>o«(0 , 

(42) r e 2 (t, x, £, s, y, r,) = -(V, + u(?7)$)y e (t -s,x-y) 

- v{() A (ufa) A V x Y t {t - s ,x-y)), 

m 2 e (t, x, f , y, 77) = -<9 t V x Y e * a G(t, x - y) . 

Therefore the characteristic flow Z e := (X e ,E e ) of the regularized 
Vlasov-Maxwell system (17) satisfies a mean- field integro-differential 
system of the form (32), where the integral operator K e is defined by 

(41) -(42). That the integral kernel of K e satisfies (30) is obvious since 
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Y e (t,x) = for all x G R 3 if t < 0. That it also satisfies the Lipschitz 
condition (31) is equally obvious in view of (26). Henceforth, we denote 

(43) L(T, e) = Lip T (/c e ) 

where k € is the integral kernel of the operator K e defined above. 

Bringing together the results obtained in the previous sections, we 
finally obtain the regularized Vlasov-Maxwell system (21) — or equiv- 
alently (17) as the mean field limit of the iV-particle system (27) as 
N ->■ +oo. 

Theorem 5.1. Let e > and f in G L°°(R 3 x R 3 ) be such that 
f m > a.e. on R 3 x R 3 , supp(/ m ) is compact, 

and 

[[ r(x,t)dxdt = i. 

(In other words, f m is an essentially bounded probability distribution 
on R 3 x R 3 .; 

Let (f e ,u e ) be the solution of (17) with initial data (18) — or, equiv- 
alently (f e ,E e ,B e ) the solution of (21) with initial data (22). 

On the other hand, let (xfjv, Cih)i<i<N be a sequence (indexed by 
N) of 2N -tuples in R 3 such that the initial empirical measure of the 
N -particle system 

1 N 

f N := j^Yl 5 *Zn ® 6 ^m f mdxd t weakl v 

i=l 

as N —>• +oo. Then 

a) for each N > 1, the system (27) has a unique solution defined for 
allt >0 and denoted t h-> (xi t N(t),£i,N(t))i<i<N ■ 

Denote 

1 N 

i=l 

the N -particle empirical measure at time t for each t > 0. Then 

b) for each T > 0, 

f N , e (t, •) -> f e (t,x,£)dxd£ weakly in ViB? x R 3 ) 
uniformly in t G [0, T] as N — >■ +oo; more precisely 

SUp distM^R (fe(t, ; ■), f N ,e(t, ; •)) < C(T, e) distMi^ (f", ft?) 

te[o,T] 

where C(T,e) = (1 + TL(T, e))e T2L ^ with L(T, e) as in (43). 



26 



F. GOLSE 



Proof. As we have explained before stating Theorem 5.1, the regular- 
ized Vlasov-Maxwell system (17) with initial data (18), or equivalently 
(21) with initial data (22) is a special case of (32) with integral operator 
K e defined by (41)-(42). 

By Proposition 4.1, the system (32) with integral operator K e has a 
unique solution 

(X e ,E e ) = (X e ,E e )(t,x,C;f in ) 
for any initial Borel probability measure f m on R 3 x R 3 . 

On the other hand, if the map t h- > £i,jv(£))i<i<jv is C 1 on 

R + , a straightforward computation shows that it is a solution of (27) 
with initial data 

(£i,Jv(0),£i,Jv(0))i<i<JV = (x™Ni Q™N)l<i<N 

if and only if 

Xi,N\t) — X € (t, X™ N , Qjf] f™) , £i,N\t) = ^e\tj X ^N^ £i™Ni Pn) ■ 

This is in turn equivalent to the fact that the empirical measure 

1 N 

/iV,e(V,-) = M^txiMt) ® S t 



1=1 

satisfies 



fN,e(t, •, •) = (X € , S £ )(i, •, •; /™)#/™ ■ 
Therefore the existence and uniqueness result in Proposition 4. 1 trans- 
lates into statement a) in the theorem. 

Comparing this last equality with (40) and using the Dobrushin in- 
equality in Proposition 4.3, we arrive at b). (We recall that the dis- 
tance distMKR metrizes the topology of weak convergence on the set 
"P(R 3 x R 3 ) of probability measures on R 3 x R 3 .) □ 

Theorem 5.1 can be viewed as a quantitative statement about the 
continuity in the weak topology of probability measures of the solution 
of the regularized Vlasov-Maxwell system (21) in terms of its initial 
data. Its corollary stated below explains how this information can 
be translated into uniform convergence of the iV-particle regularized 
electromagnetic field to the one created by the solution of (21). 

Let (E e , B e ) be the electromagnetic field associated to the solution 
of (21) with initial data f m and let (E N , e , B N>e ) be the solution of the 
regularized Maxwell system 

div x B Nj€ = , 
d t B Nt£ + curl,, E Nte = , 

dlV x E N>e = Xe *x Xe *x Pf Nit , 

d t E Njt - curl x B Njt = -Xe *x Xe * x 3f N , e , 



(44) 



MEAN-FIELD LIMIT FOR VLASOV-MAXWELL 27 

with initial data 

(45) E N,e\ t=0 = -VxXe*xXe*xG* x Pfin , B e \ t=Q = , 

In other words, (E Nte , B N>e ) is the regularized electromagnetic field cre- 
ated by the iV-particle system. 

Corollary 5.2. Under the same assumptions and with the same nota- 
tions as in Theorem 5.1 
a) for each T > 0, one has 

dist M xjj {Pfe(t,-),p fNit (t,-)) < C(T,e)dist 

MKR (/' j /jv ) 

and 

~ W*>-)lk-M ( R3) < 4C(T, e) disW {r,f%) 
for all t G [0, T] and N > 1 — so that in particular 

1 N 1 N 

xJ2 5 *iMt) -> Pfettr) and — -> j/ e (t, •) 

i=l i=l 

uniformly in t G [0,T] as N ^ +oo, where the first convergence holds 
in the weak topology ofV(R 3 ), while the second holds in the weak topol- 
ogy of bounded Radon measures on R 3 / 

6j /or eac/i T > 0, (-Ejv,e, 5jv )t ) converges to (E e ,B e ) uniformly in 
(t, x) G [0, T] x R 3 with the following estimates 

\\B e - -BAr, e ||L°°([o,T]xR 3 ) < C"(^, e) distM/s:R(/ m , /™) 

and 

ll-E'e — -£^v,e||L°°([o,T]xR 3 ) < C"(T, e) dist A f^ R (/ m , f^) , 

w;/jere 

C'(T,e)=4TC(T, e)\\W x Y £ \\ 

and 

+\\d t V x Y e * x G\\ 
and where C(T,e) is as in Theorem 5.1. 

Proof. For the first statement in a), pick any bounded and Lipschitz- 
continuous function on R 3 and observe that 

/ (Pf e (t^) - p fN}C {t,z))<f){z)dx 

</R 3 

(x) - </>(y))n(t, dxd^dydn) 
< 2\\<b\\ ^/i,oo / / 1 A \x — y\7r(t, dxd^dydrj) 

</</R 3 xR 3 xR 3 xR 3 



R 3 xR 3 
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where n(t, •) is any probability measure in U(f e (t, •), fN,e(t, •))■ Mini- 
mizing the r.h.s. in 7r(£, •) while maximizing the l.h.s. in such that 
Lip(0) < 1 and applying the Kantorovich duality theorem (see formula 
(5.11) on p. 60 in [19], or Theorem 1.3 on p. 19 in [18] and Remark 
7.5 (i) on p. 207 in that same reference) leads to 

<Hst M KR(Pf e (t, ■), PfN^i ■)) < distMKR(fe(t,-)jN,e(tr)), 

and one concludes with Theorem 5.1 b). 

For the second statement, proceed in the same way, replacing the 
test function with if)(x) ■ v(£), assuming that ip G W^ 1,00 (R 3 ; R 3 ). 
Thus 

/ UfAt, z) - j fNit (t, z)) ■ ^(z)dx 
= / / (^0*0 • u(0 - i>(y) ■ v(rj))ir(t, dxd^dydr]) 

J </r 3 xR 3 

< (2||0|| w i,oc + 2||V||loc) 
x (1 A \x - y\ + 1 A |£ - r)\)ir(t, dxd^dydi]) . 

J ./R 3 xR 3 xR 3 xR 3 



Minimizing the r.h.s. in ir(t, •) while maximizing the l.h.s. in ip such 
that II^IIvk 1 ' 00 ^ 3 ) < 1 leads to 

~ j/ JV , e (*,-)lllV-M(R3) < Ad38t MKR (f e (t,-),f Nte (t,-)) 

and one concludes again with Theorem 5.1 b). 
Now for b). Observe that 

B e (t,x)-B Nie (t,x)= / V x Y e (t-s,x-y)A(j fe (s,y)-j fNe (s,y))dy 

JO JR 3 

so that, for all x G R 3 , one has 

\B e (t,x) - B N:t (t,x)\ <t sup ||V x y e (s,-)||wi.<»(R3) 

0<s<t 

x sup \\j fc (t,-) - Jf N .Mr)\\w-^(^) 

0<s<t 
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and one concludes by a). The estimate for the electric field is obtained 
in a similar way, from the formula 

E £ (t,x) - E Nj€ (t,x) 

= - [ [ d t Y e (t-s,x-y)(j fe (s,y)-jf Ne (s,y))dy 

Jo JR 3 

- / / V x Y e (t- s,x -y)(p fe (s,y) - p fNe (s,y))dy 
Jo Jr. 3 

- / (V x d t Y £ (t, ■)* x G)(x-y)(p fi n(y) - Pfi Jy))dy. 

</R3 

□ 



6. Final remarks and conclusion 

The discussion above leaves aside several additional questions about 
this regularized dynamics. 

First, although we insisted on choosing Rein's regularization proce- 
dure because the "pseudo-energy" 

W e ](*):= // e(Z)f e (t,x,Z)dxdt+ f l(\E e \ 2 + \B e \ 2 )(t,x)dx 

</</R3xR3 JR3 

is constant under the dynamics of (21) — see Proposition 3.1 — we 
have not used this quantity in the mean-field limit itself. (The same 
can be said of earlier work in the same direction, first and foremost [3] 
and [7].) 

Strictly speaking, Rein's theorem in [16] establishes that W[/ e ] is 
constant only in the case where the initial particle distribution func- 
tion f m is a probability density (in fact, in the more general case where 
f m > a.e. is a measurable function such that W[/ m ] < +oo.) There- 
fore, even though the empirical measure /jv, e in Theorem 5.1 is a weak 
solution of the regularized Vlasov-Maxwell system in the sense of mea- 
sures, Rein's theorem cannot be applied directly to initial data of the 
form /jy . 

Proposition 6.1. Let N > 1 and e > be fixed, let (x- n , £- n )i<i<Ar be 
a N -tuple of elements in R 3 x R 3 ; and let 

1 N 

i=l 
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Let t i — y (xi(t),£i(t))i<i<N be the solution of (27) with initial data 
(^r^i n )i<i<JV; and define 

1 N 

fN,e(t, ; •) := — <W) ® ' t ^°- 

i=l 

w[/^(t) = wra 

/or oZZf > 0. 

Proof. Regularize /jy in (x,£), for instance by replacing f™ with its 
convolution f™ * X £ ( v with any (nonnegative, compactly supported) 

mollifying sequence (C;)r/>o on R 3 x R 3 - If \%T\-> 10" I < ^ f° r eacn 
i — \,...,N and supp(C) C -8(0, rf) x 5(0, 77), one has 

supp(/£„) C B(0, i2 + V ) x B(0, R + V ). 

Then 

E™ ;e = —V X G * x Xe *x Pfi* and E™^ tt = —V X G * x Xe * x Pffr 

both belong to L 2 (R 3 ) since V x G(x) = 0(\x\~ 2 ) as \x\ — > +00. 

Denote by (/jv j7?je , £jv,??,e> B^ , q ^) the solution of (21) with initial data 
{f l N,rf> Xe*x ^ l N,t],ei 0)- Applying Rein's theorem (Proposition 3.1) shows 
that 

W^, e ](0 = W^]<+oo, t>o. 

Let 77 — >■ 0. Then fff — > f™ in the weak topology of Borel proba- 
bility measures on R 3 , and 

/ / e(t)f$ tV (x, i)dxdi ^ ff e(0f N n (dxd0 

JjR 3 xR 3 iJR 3 xR 3 

since supp(/#^) C B(0,R + rj) x B(0,R + 77). On the other hand, 
denoting 

which is a regularizing sequence on R 3 , we see that 

rpin n fjiin i Z?* ra 

N,T),e — U r) *x &N,e N,e 

in L 2 (R 3 ). Therefore 

w#„] mm < +°o as v o . 

Reasoning as in the proof of Theorem 5.1 and Corollary 5.1, since 
fN, v ,e and /jv, e are two solutions of the same regularized Vlasov-Maxwell 
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system with initial data satisfying f]tf — > f™ in weak-"P(R 3 x R 3 ) as 
77 — > 0, we deduce from Dobrushin's estimate that 

f NtV , e (t, ■, •) ->• /jv, e (t, •, •) in weak-P(R 3 x R 3 ) 

uniformly in t G [0, T] for each T > 0, while 

uniformly on [0,T] x R 3 for each T > as 77 — > 0. (Recall that 
(E NjJ):e , B NiV>£ ) is the solution of the regularized Maxwell system (19) 
with initial data (20) with fN,r),e m the place of f e , while (E Nj€ , B Nje ) is 
the solution of the same system with f^, e m the place of f e .) 

Obviously supp(/jv j7? , e (i, •, •)) C -8(0, i? + 77 + £) x R 3 (since particles 
travel at speed < 1). Since 

/ PfN :V A^ x ) dx = 1 and / \jf N ^ e (t,x)\dx < 1 

we conclude from the formulas expressing E N>v ,e and B N ^ e , i.e. 

^jv,jj,£ = — <9tV x F e * X G* X Pfin — V x Y e * t)X Pf Nr) e — d t Y e * tjX jf Nri e 

B N>Vie = cml x {Y e * tjX j fNve ) 

that there exists C 6: t > such that 

\E NjV , e (t,x) \ + \B N ^ e (t,x)\ < C 6:T 

uniformly in 77 > as (t, x) G [0, T] x R 3 . Therefore, for each t G [0, T], 
one has 

supp(/^, e (*, •, •)) C 5(0, R + rj + T) x 5(0, i? + 77 + 2TC„t) 
uniformly in 77 > 0, and 

suppCp/^C*,-)), supp(j /jVi)e (t,-)) C 5(0, i? + 77 + . 

We also deduce from the formulas expressing .EW^e and B^^e that, if 
< t < T, then 

-Eiv,?7,e — <9 t F * x E™^ e and -Bjv^e 

have compact support in [0, T] x 5(0, R + 77 + e + 2T). 

Since we already know that 5)-?^ — > E™ in L 2 (R 3 ) and that 
(E N ^ e ,B N ^ e ) ->■ (E Nj£ ,B N:e ) uniformly on [0,T] x R 3 for each T > 
as 77 — > 0, we conclude that 

W[fN, V ,e](t)^W[fN,e](t) 

for each t > as 77 — > 0, so that 

mfN,.](t) = wir N n j 

for all t > 0. □ 
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Since the quantity W[fN,e](t) is an invariant of the motion for the 
dynamics of (27), a natural question would be to check whether (27) 
is a hamiltonian system. We have left this question unanswered for 
lack of applications of immediate interest. Besides, since the molli- 
fied dynamics (27) is not a fundamental law of physics, whether it is 
hamiltonian is perhaps a purely academic question. 

Several remarks are in order to conclude the analysis presented above. 

First, we have chosen to study (27) instead of (28) because empirical 
measures constructed on the trajectories of (27) are exact solutions 
of (21) in the measure sense. This interpretation disappears if one 
neglects the self-interaction term in the regularized system. However, 
because of Proposition 3.2, both regularized systems are equivalent for 
e > fixed and N — > +oo. On the contrary, none of these systems 
seem to make very much sense for e = 0, and in any case modeling 
self-interaction in classical electrodynamics is a source of conceptual 
difficulties, as mentioned above. As explained on p. 677 in [8], there is 
little hope of deriving the Vlasov- Maxwell system as a mean-field limit 
of a particle system without regularization as long as self-interaction is 
not better understood. 

Also, the work presented here leaves aside several issues mentioned 
in [8] without further analysis, such as including spin in the particle 
system; as in [8], the regularization procedure used here (based on 
convolution with a blob function in the space variable only) is not 
Lorentz invariant. In any case, this regularization should not be viewed 
as an attempt to depart from a strict point particle model and introduce 
some physical effect involving particle size, but as a mere mathematical 
expedient. 

As a matter of fact, the regularized dynamics (27) could be regarded 
as an analogue for the Vlasov-Maxwell system of the well-known vor- 
tex blob method used in numerical simulations of incompressible fluid 
flows: see section 5.3 in [14] for a concise, yet lucid presentation of this 
method. In the present context, one has the following approximation 
result. 

Proposition 6.2. Let f m G L°°(R 3 x R 3 ) be a compactly supported 
probability density on R 3 x R 3 . 



For each N > 1, let {xf N ,£™ N )i<i< N G (R 3 x R 3 )^ be such that 




i=i 



in the limit as N — >■ +oo. 
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For each e > 0, letter (x(t)i t N,e,^(t)i,N,e)i<i<N be the solution of the 
regularized dynamics (27). Then, there exists subsequences N k — > +oo 
and e k — > such that 

1 N 

jyE^,^^W 6 ^N k ,, k (t) ->• /(*, 0<*cd£ ^ea% m P(R 3 x R 3 ) 

i=l 

uniformly int<E [0, T] /or eac/i T > 0, where (/, E 1 , 5) is a global weak 
solution of the Vlasov- Maxwell system (1) with initial data (2)-(3). 

The global existence of weak solutions of the Vlasov-Maxwell system 
was obtained by DiPerna-Lions [6], with some additional precisions to 
be found in [16]. The (very weak) approximation result above follows 
from applying Rein's Proposition 4 in [16] and Theorem 5.1. Indeed, 
let f e be the solution of the regularized Vlasov-Maxwell system (17) 
with initial data (18) for each e > 0. 

By Proposition 4 in [16], there exists a subsequence e k — > such 
that fe k (t, •, -)dxd£, — > f(t,-,-)dxd$, in weak-P(R 3 x R 3 ) uniformly in 
t G [0, T] for all T > as e k — > 0. Therefore, for each k > and T > 0, 
there exists e k > small enough so that 

tist MKR (f(t, •), fe k (t, •)) < 2~ fc - 1 for each t E [0, T] . 

Then, for this value of e k > 0, by Theorem 5.1, there exists N k > 
large enough so that 

dist MKR (f ek (t, - , •), /iv fciefc (t, •, •)) < S'"- 1 for each t G [0,T] , 
where 

1 N 

fN k ,e k (t, ; •) := ^2 S xi,N k ,e k (t) ® %,iv fc , £fc W • 
i=i 

One then concludes by the triangle inequality. 

Finally, it could be interesting to study fluctuations as Braun and 
Hepp [3] did for the Vlasov-Poisson case; this is left for future investi- 
gation. 
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